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SUMMARY 

A study has been made of certain features of transverse resonance 
in an idealized combustion chamber. Modes of weak oscillation are des- 
cribed^ with avereige Mach number neglected hut ^^ith the effect of axial 
gradient of mean temperature considered. These weak oscillations are 
assumed to he aaqplified by coupling with the vigorous combustion occur- 
ring in flame-holder wakes. 

A flame -holder wake is idealized as a thin annul ns performing 
mechanical work on the surrounding gases. The amplification of the 
various modes depends on the flame -holder diameter, centerhody diameter, 
and on time-lag effects. 

The nature of finite transverse periodic waves is also analyzed. 
Results show that such mves have frequencies independent of angplitude 
end do not steepen with time, at least to second order. 


INTRODUCTION 

The flow in the combustion ctsmhers of jet engines is frequently 
subject to a reso n a n t oscillation, the character of which depends on the 
paxticular engine configuration. The transverse -wave motions described 
by the classical theory of acoustic oscillations inside cylindrical 
enclosiires (refs. 1 and 2) may sometimes be inqportant. 

In the classical theory, the fluid medium is considered to be at 
rest and in a uniform state except for weak isoenergetic fluctuations 
in velocity and state properties. Under these assumptions, the motion 
is governed by the wave equation, written here in cylindrical coordinates 
for the pressure disturbance (sketch (a)) 




+ p' 
j^rr 
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Velocity and pressure fluctmtions are related by the equation 

P + 7P’ = 0 (lb) 

A class of solutions of equation (la) suitable for internal flows is 



where n is any positive integer. 



(a) Coordinate system and sketch of fundamental sloshing mode. 

(a con^jlete list of symbols is given in the appendix.) For resonant 
oscillations, C and CD are determined by conditions at the ends of the 
cylinder and from the condition of zero radial velocity at the cylinder 
t-rall (r = R) 


J^(C) = 0 (3) 

For each value of n there exists an infinite nuniber of solutions of 
equation (3) for C, becaiise is an oscillatory function of its 

argument. Hereinafter, the number m will be used to specify that C 

is chosen as the (m + l)^^ number for which the derivative of J 
vanishes, but does not. 

If consideration is restricted to oscillations independent of x, 
the two sin^jlest cases are: first "radial" mode (n = 0, m = l), 

3.83i§t r 

p' = e ^ Jo^3.83 Ij 


(t) 
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and first "sloshing" mode (n ?= 1, m = O), 



The loci of disturhance motions for the first sloshing mode (n = 1, 
m = O) are shown in sketch (a) . There are two nodes,, diametrically 
opposed, and the fluid appears to slosh between these nodes in the illus- 
trated case of a standing wave. This mode has the lowest frequency 

— of all the possible transverse cylindrical J^aves. 

In the present report, it will be assxmied that resoneint combustion- 
chamber osc ill ations result from the amplification of initially weak 
transverse pressure disturbances interacting -v/lth the combustion pro- 
cess . (The amplification of transverse waves in solid-propellant rockets 
has been studied by Grad (ref. 3) and Chdng (ref. 4).) On this assimption, 
the follovring problems are investigated: 

(1) Modes of osc il lation; The transverse modes of oscillation shoiild 
be related to conditions in a combustion chamber, which contains gases in 
motion with severe gradients of state properties o'vring to heat released 
by the conbtistion process. These conditions are at variance with the 
assumptions of classical acoxistic theory. 

(2) Amplification at a thin annulus; Flame holders (fig. l) induce 
especially vigorous burning in their wakes. Flame-holder configurations 
would therefore be expected to affect the amplification of combustor 
oscillations in the various modes. Amplification is also affected by any 
time lag in the interaction of the burning process with the oscillations. 


( 3 ) Finite transverse waves; Combustion-chamber oscillations may 
have extremely high amplitudes. The character of finite-strength trans- 
verse waves therefore merits an Investigation to detenmine, for example, 
whether or not compressive transverse waves steepen to form shocks, as 
do plane waves. 


MODES OF OSCIUATiaN 
Differential Equations and Assumptions 

The action of a flamie holder in a combustor (fig. l) tends to pro- 
duce a region extending downstream in which burning is especially vig- 
oro'us. The gases surrotmding the flame-holder wake would presumably be 
bvuming, but at a slower rate- These surrounding gases (shown shaded 
in fig. 1 ) are assumed to obey the equations of compressible nonviscous 
flow -vrLth heat addition 
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Momentum: 

Continuity: 

Energy: 

State: 


Dg 




+ V . (pq) = 0 


Dc^ 

P = (r - i)pc^ 



( 6 ) 


In tte atsence of distirrbance, one-dimensional steady flow is 
ass\jmed, and equations (6) become 


= - Px 
pu = constant 

p = (r - l)pc^ ^ 


( 7 ) 


The disturbance field is assumed to involve the high-frequency 
motion \xsual in cases of acoustic resonance. The consequent assumption 

^ = order^S ^ (8a) 

is made concerning derivatives in the distTorbance field. Combiostors 
usually produce large temperatTire changes at rather low Mach number. 
Accordingly^ the eissim^ition of small Mach number is made 

M « 1 (Sb) 


The various flow quantities may now be assumed to be pertiarbed 
slightly from their steady values; thiis, for exaji5)le, 

p(x,r,0,t) = p(x) + p'(x,r^e,t) (8c) 

If eqmtion (8c) and similar expressions for the other flow prop- 
erties are substituted into eq'uation (6), taking note of equations (v)^ 
(8a), and (8b), there follows, to first order. 



^ p' - v2p> 
T -^x 
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If it is further assumed that, when the temperature undergoes a relative 
change, the rate of heat release H has a relative change of the ^ame 
order, then the right side of the foregoing equation is of order M, and, 
for consistency with equation (Sb), should he neglected. The equation 
replacing the mve equation (la) is thus 




(9a) 


To the same order, the following additional equations apply: 


(9b) 



(9c) 


^ ^ ^ (9d) 

p p T 

Thus, to zero order in Mach nimiber, the osci ll ations are not ampli- 
fied by coupling between an oscillation and the coirbustion process. How- 
ever, to this order in Mach mmiberjj_ the mode of oscillation is affected 
by the mean temperature gradient Tj^ produced by the burning. 

Solutions 

Equation (9a) may be satisfied by solutions in the following form 
(obtained by separation of variables); 




i) + ^>'n(c i) 


where G(x) must satisfy the equation 


G"(x) + ~ G'(x) + 
T 


CD 


a^(x) R*' 


G(x) = 0 


(lOa) 


(10b) 


If the distance reqtiired for \mit relative change in temperature is long 
conpared 'vrLth the bu^er diameter, then the term of equation (lOb) con- 
taining the factor T^ may be neglected, and equations (lO) may be re- 
placed by the classic^ solution (eq. (2)), the speed of soimd 
a = -y(y-l)cpT being regarded as variable with x. If the term involv- 
ing T is not neglected, then equation (lOb) must usually be integrated 
numerically. 
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The constant C is determined hy transverse "boundary conditions ^ 
and CD is then determined hy an axial boundary condition, 'which is 
assumed for this case to be 


CD 



(lOc) 


■where a^ is the final (maximum) speed of sound attained in the com- 
bustor. If equation (lOc) applies, the x-dependence of the oscillation 
is governed by the eqiiation (which replaces eq. (lOb)) 


G" 



T 





G = 0 


(lOd) 


Clearly, if T approaches a constant value as x increases, eqtiatlon 
(lOd) provides that G must become and remain constant (G "being 
assxmed to remain finite). Therefore, the effect of equation (lOc) is 
to assume that any osci ll ation becomes entirely transverse for large 
X. In the region ;q)stream of the flame holder there is a sinusoidal 
variation of pressure 'vri.th x, ^diich for a traveling ■wave results in a 
spiral ■^•ra.ve form. 


Figure 2 illustrates the axiaJ. form of the mode corresponding to 
equation (lOc) obtained by nimerical integration of equation (lOd) for 
the illustrated axial veiriation of sound speed and for the ratios 
R/L = 0.1 and 1. The constant C 'tiias taken to correspond to the first 
sloshing mode (n = 1, m = O), and the solutions were adjusted so that 
G(«») = 1 in each case. The larger -value of R/L corresponds most 
closely to a large combustor, while the other -value would be more reason- 
able for a combustor with a very small diameter. 


AMPLIFICATION AT A THIN ANNULUS 
Assun^jtions and Boundary Conditions 

It is assumed that burning in the -t-Tale behind a flame holder, 
here assumed to be a nnu lar, is substantially more -vigorous than in the 
surrounding gases (fig. l) . Thus, the assim 5 )tions leading to equations 
(10) are inappropriate to oscilla-tions occurring -within the annulus. A 
different sin 5 )lifying assumption is made concerning the annulus j namely, 
that the region of vigorous burning is extremely thin - too thin to 
support a transverse pressure difference. Thus, en^ploying subscripts c 
and w to denote the dis-turbance fields between centerbody and annulus, 
and annulus and outer wal.l, respectively, 

PcC^l^e^t) = p^^(r^,e,t) 


(Ha) 
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Within the a nn ulus, the burning rate is assumed to undergo relative 
change in proportion to relative change in temperature (and hence pres- 
sure, for sufficiently small axial steady-state gradients (eq. (9)). 

This fluctuating energy release in the annulus is communicated to the 
surrounding gases as fluctuating mechanical vork, such as trould be per- 
fomned by an aneroid expanding and contracting. !Ehe mechanical work 
would be in proportion to the radial velocity difference 

- v^(r^,0,t) 

Thus, in view of equation (9b) and the postulated relation betwen tem- 
perature and heat release, for a vanishingly thin annulus, 

^ ' 5p ' ^ ' 

^ (r^,0,t) - (rj_,e,t) (lib) 

where K is a qiiantity -v/hich depends on the combustion process tjlthin 
the annulus, and must In general vary ^Tith x. If the annulus bnp an 
axial extent much greater than the combustor diameter, then K may 
be regarded as constant. (This condition is not likely to be satisfied 
in a real combustor, but is considered to be qualitatively permissible.) 

Equations (Ua) and (Ub) (K regarded as constant) are boundary 
conditions to be inposed on the oscillation in the surrounding gases. 

This oscillation is considered as composed of two parts, applicable for 
r < r^_ and r > rj_ separately, each part having the form of equation 

(lOa) . Thus, two sets of constants appear; etc, and o^., Further 

boundary conditions specify vanishing normal yelocity at the outer wall 
and at the centerbody: 

^ (R,e,t) = 0 (He) 

^ (ro,0,t) = 0 (Hd) 

The foin* homogeneous boundary conditions (eqs. (ll)) relate the 
four constants a^, p^, o^, and p^ and yield the following determinantal 
equation for O): 
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Ji(Z) 


^There z = 2r/a^^ Z = 2R/a^, and ZQ = ffirQ/a^. First-order terms in K 
define i] and yield the amplification, which is eq\ial to the negative 
of the imaginary part of m: 


Im(aj) 

jcr3_Ka^ 



(13) 


The grouping of constants on the left should he held constant for com- 
parison of the effect of changes in flame -holder and centerhody config- 
uration, In particxfLar, since the quantity K (eq. (lib)) must depend 
on the rate of heat release per unit surface area of the annulus, rtr^^K 
is proportional to the total rate of heat release at the annulus per unit 
length, nrit^ hence, for comparison purposes, should be held constant. 

If there is no centerhody. 


Im(cD) = z2 

jtrj_Ka^ 






( 14 ) 
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The amplification k^, for the case of no centerhody, is plotted in 
figure 3 for several inodes as a function of r/R. The greatest an 5 )lifi- 
cation occurs i-ri.th the radial mode (n = 0, m = l) -v/hen the burning is 
concentrated near the axis. Angplification is generally greatest when 
energy is added near a pressure loop (ref. Z, p. 226), and only the radial 
modes (n = O) have such a loop near the axis. The higher modes (n > l) 
have pressure nodes at the axis and are not an^lified at n.1 1 when burn- 
ing is concentrated there. If the annulios is placed near the wall, higher 
modes (which have pressure loops at the wall) undergo the greater 
an^illf ication. 

The effect of a centerbody is shown in figure 4. Only the first 
sloshing mode (n = 1, m = O) is considered. If the a.Tvrmlus is near the 
outer ■viall, the presence of a centerbody decreases the amplification 
somewhat, I'diereas the opposite is true if the annulus is near the center 
of the chamber. 


& 


Time Lag 


Figures 3 and 4 predict the occurrence of the radial mode for burn- 
ing near the axis or higher-order transverse modes for burning far from 
the axis. However, it is reasonable to inquire whether or not these re- 
sults might be altered if there is introduced a time lag betwen state 
changes and fluctuations of the combustion process. The slmrpiest lag 
hypothesis for the present problem is that a change in temperature pro- 
duces a corresponding change in heat release at a somewhat later time. 
(Crocco and Cheng (ref. 5) and Cheng (ref. 4) consider the effects of 
time lag in their analyses of resonance in rocket engines . ) The basic 
cause of such a lag is not specified, but might derive from the ccmbus- 
tion kinetics of each burning element. 


Thus, equation (Ub) is replaced by 




3p ' 

^ 


(r^, 0,t) = K (r^^, 0,t-T) 


(15) 


where t is a (constant) time lag. The analysis subsequent to equation 
(Ub) may then be repeated, resulting in a modification of equations (l3) 
and (14) 


5 Im(cD) = [A, or A_] cos S2r 

:tr^Ka| ^ 


( 16 ) 
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If 2 tt < Increasing the frequency decreases the amplification and^ in 
fact, causes damping if it/2 < fir < it. Therefore, the an 5 >lification of 
the first sloshing mode (n = 1, m = O) trould he more intense relative to 
the other modes than is indicated in figures 3 and 4. 


Actually, it tTOuld not he expected that a single time lag t char- 
acterizes the entire huming process, hut rather that each- burning ele- 
ment trould have a different response lag, so that each of the quantities 
A cos 2 T in equation (l6) is replaced by 



cos 2 tt dT 


■vThere 

be 


F(t) is a distribution function for which a suitable choice might 


F(t) 



(17) 


which has the properties 

has a maximum at t = tt* 
place of equation (l6). 



Use of the preceding 


F(r)dT = 1, and which 
relations yields, in 


— Im(a)) = [A^ or A^] 
jcr^^Kaj 


1 - (2T*)^ 


( 18 ) 


Therefore, modes for which 2 t > l n-TOuld be danced if this equation 
holds. For an^Jlifled modes, the coefficient of the A's diminishes the 
an^ilification of modes of higher frequency. 


FIDJITE TRANSVERSE WAVES 

In one -dimensional unsteady flovr, strong con 5 )ression waves progres- 
sively steepen into shocks, •(•?hich move ’(■rith a speed greater than that of 
sound. This dependence of propagation speed on an^litude ingplies a de- 
pendence of \Ta-ve form and frequency on amplitude in any one -dimensional 
cyclic process involving compression \raveB. The subsequent analysis is 
concerned with determining whether or not amplitude is similarly impor- 
tant for strong transverse oscillations in a cylinder, in order to indi- 
cate whether agreement between acoiistic theory and experiment may be 
anticipated. 
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Differential Eqmtions 

An isentropic irrotational disturbance field (not necessarily -weak) 
obeys the equations 

.2 .. . -2 ^ 


Momentum: 


«P 4 - +V + 


jr £ a_ 


*'t "2"^r-ip“r-i 

Continuity: + 7-(p£) = 0 

State isentropy: p <x p^ 

Definition of potential: q = y cp 


(19) 


J 

■v^ere a is the speed of sound in the absence of any oscillation. Equa- 
tions (l9) may be combined to yield a differential equation for cp : 


—2 2 
(Ptt - a V <P = 


<p = - + \ a-v(<l^) + (r - l) f (20) 


Next, a finite-strength vrave pattern is sou^t in the form of a 
power series in an amplitude parameter e 


cp = aR 


(ecp^^^ + e^cp(2) ^ 


( 21 ) 


Substitution into equation (20 ) and collection of terms of first and sec- 
ond order in e yield, respectively, 

- a^V^cp^^) = 0 (22) 


and so forth. 


Solutions 

One -dimensional motion . - A weak traveling wave may have the foim 



= k sin S2 



(24) 
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•whicti satisfies equation (22). Then equation (23) may he -written 

_ (2) -2 (2) r + 12% ,2 . . x\ 

'Ptt - ^ 'Pxx - — T" 22^t + =j 

Ko solution of equation (25) can he -written consistent with the existence 
of a strong -wave moving at constant speed -with a permanent wave form, 
heca-use the solution must contain a particular integral of the type 

- Y ^ ^ + =) 

Therefore, no permanent one-dimensional finite -wave is possible (ref. 2, 
p. 480). 



Transverse motion in a cylindrical pipe . - From equation (2), the 
first sloshing mode (n = 1, m = o) is chosen for study: 

(p(l) = J^(z) sin (st + 0) (27) 

Eq-uation (27) represents a wave spinning in a clockwise direction. The 
addition of a -^rave spinning in the opposite direction -wo-uld yield a 
standing, or sloshing, motion. Equation (27) cannot he interpreted as 
representing a succession of piilses traveling with the speed of sound, 
as can eq-uation (24). The "wave fron-t" of equation (27) travels at an 
apparent Mach nuiriber -varying from zero at -the cylinder axis to 1.84 at 
the cylinder wall and therefore should he regarded essentially an an 
interference pattern rather than as a -wave front. The foregoing consid- 
eration suggests that the progressive steepening -tvith time, characteris- 
tic of isolated pulses and ill \istrated hy equation (26), need not neces- 
sarily he expected in the case of treinsverse motion in a cylinder. The 
attenipt is therefore made to obtain in the foim 

cp(2) = 

sin 2(S?t + 0) (28) 

Substituting equations (27) and (28) into equation (23) yields an equa- 
tion for f (z) 


f.. + ^ + 4^1 - ^ 


(29) 
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The houndaxy condition at the cylinder wall is 

f'(z) = 0 

and^ at z = 0^ f(z) must he an analytic function. 
The solution to equation (29) may be written as 


(30) 


f(z) = I |(r + 1 ) [dJ2(2z) - g(z)J + j 

where^ by substitution into equation (29), g(z) is a partic\ilar solution 
of the eqmtlon 

g" + i g' + 4^1 - (32) 

and DJg(2z) is the con5)lementary solution of the same equation. The 
constant D is chosen so that boundary condition (30) is satisfied. 
Eqxjation (32) may be solved by a series method, adopted for convenience 
in computation, yielding the result 



7 

187 /zN^® 437 fz\^ 

24\2j ^ 180^2 ) 

8640 \2 ) 151200 \^2 j ~ 

1979 1 

rz_Y^ 4387 / 


7257600' 

(2 j ^ 228614400 \ 

¥) - ••• 


(33) 


Introduction of equations (3l) and (33) into boundary condition (30) 
gives 


D = - 0.384433 (34) 

Equations (28), (3l), (33), and (34) complete the solution, indicating 
that a "permanent" transverse an oust ic wave of finite amplitude is pos- 
sible, at least to second order. Thiis, a strong wave form may have a 
Fourier series, consisting of the linearized solution and its harmonics, 
for transverse motion in a cylinder, but not for a purely one-dimensional 
\T&ve. 


Numerical Results and DiscTission 

The resTilts of the previous section may be expressed throu^ the 
equations of motion (eqs. (l9)) in terms of pressure fluctuation: 
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=0 


1 + 6*Pj^ cos 


where hhe definition 


+ 0) + e^^jP2o + 1*22 2(SJt + 9^ 


+ • 


® ■ rzjj^iz) 

is made^ so that Pj^(z) =1, an a matter of choice. Furthennore, 



With y = 1.4 assumed^ numerical results for ^20.> ^22 

are shown iu figure 5 as functions of r/R (equivalent to z/z). The sym- 
bol Pj_ represents the radial distribution of pressure anplitude in a 
spinning acoustic vmve of the first sloshing mode (n = 1, m = O); PgQ 
is a steady second-order correction to the radial pressure distribution^ 
and signifies that the oscillation tends to exhaust the center of the 
cylinder and congjress the fluid near the wall, in the manner of a cen- 
trifugej and P 22 is the radial distribution of pressure amplitude of 
the first overtone. This overtone is strongest relative to the funda- 
menteil tone at the wall and is relatively ■(•reahest at the cylinder axis. 


The modification of the wave form, when the oscillation is rather 
weak but not of infinitesimal strength, is illxistrated in figure 6. 

Ovring to terms of second order in e* the pressure crests at the cyl- 
inder VTall are steeper and higher than those of the fimdamental harmonic 
wve, and the trouts ane broader and shallower, in the familiar pattern 
of second-harmonic distortion. A lesser second-order effect is the 
slight increase in average pressure at the cylinder wall, owing to F 2 O' 

The steepness of the crests illustrated in figure 6 is suggestive 
of the steepening of one-dimensional compression \raves leading to shock 
formation, but is not a progressive steepening -^rLth time as in the one- 
dimensional case. Rather, the modified wave sho>m in figure 6 may be 
imagined to proceed around the cylinder -tTall without change ad infinitum. 
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To the present order of approximation, the fundamental frequency does 
not change either with time or with an 5 )litude, although overtones appear 
in proportion to arrplitude . It may therefore be concluded that violent 
comb;istion -chamber resonance is not necessarily incon^iatible with 
acoustic -theory predictions. 

The permanent isentropic character of the solution obtained in the 
foregoing analysis suggests that strong transverse oscillations may not 
be subject to the dissipative normal stresses which tend to limit the 
an 5 )litude of a one -dimensional strong-wave cycle involving shock waves. 


CQNCnJDDfG EEMARKB 

On the assumption that combustion-chamber reponance results from the 
amplification of an initially weak wave of transverse type coupled trith 
a pressure-dependent rate of energy release in the combustor, three prob- 
lems were discussed: 

1. Modes of oscillation: It was shown that if the mean flow Mach 

number is aasimned very small, then the possible modes of \reak disturbance 
are governed by a \Ta.ve eqiiatlon in classical form except for a term de- 
pending on the streamwise variation of tenperatirre in the combustor. 

Tvro such modes were determined by numerical integration, subject to the 
assumption that in the afterportion of the combustor the oscillation is 
purely transverse. These modes conform to the classical spiral modes 
(if the wave is a traveling one) upstream of the burning zone, with a 
pitch depending on the ratio of chamber diameter to length of burning 
zone, as well as on the combustor temperature. 

2. Amplification at a th in annulxis; In the approximation of very 
small Mach number, the burning process does not amplify the modes des- 
cribed. The especially vigorous burning in a flame-holder wake would 
provide amplification, hoover. A flame-holder wake is idealized as a 
thin annulus performing work upon the STirroxmding gases in proportion 
to the instantaneous pressiire level. 

AnELlysis of the res\ilting amplification indicates that concentra- 
tion of vigorous burning near the center of the combustor amplifies the 
radial mode miost strongly. Burning concentrated neax the combustor tTall 
amplifies the higher sloshing modes. The presence of a centerbody 
slightly increases (decreases) the amplification of the sloshing mode if 
energy is added near the inner (outer) wall. 

Supposition of a time lag between temperature disturbance and heat- 
release rate in the ann ulus favors amplification of lo^rer-frequency 
modes; in particular, the first sloshing mode, which has the lowest fre- 
quency of fl.l 1 possible modes . 
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3. Finite transverse waves: The nature of strong transverse acous- 

tic oscillations was considered. It was shown that, to second order, 
amplitude affects wave form, in a spinning oscillation, causing the pres- 
sure crests to he higher and sharper and the troughs shallower and 
broader. The average pressTire is dim in ished near the axis and increased 
near the wall, suggesting a sort of centrifugal pumping. The wave form 
does not undergo progressive change leading to formation of a shock, as 
in one -dimensional flowj rather, the fundamental frequency does not de- 
pend on an 5 )litude, to the present order of approximation. 

Furthermore, the permanently isentropic nature of the transverse- 
wave solution sxiggests that transverse modes in a cylindrical enclosure 
are not inhibited hy viscous effects to the same degree as are one- 
dimensional waves, and for this reason are perhaps better able to attain 
extreme amplitudes. 


Le-vjis Flight PropiiLsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 19, 1954 


3057 



mCA TN 3152 


17 


APPEHDIX - KOTATION 

The fo3JLo>ring notation is used in this report: 


a 


P 

D 

f 

G 

g 

H 

Im 

Jn 

K 

k 

L 

M 

m 

n 

^ 1 ^^ 20^^22 

P 


an5)lification rates (egs. (l3) and (l4)) 

speed of sound, - l)Cj,T 

final speed of sound attained in combustor 

constant (eqs. (2) and (lO)) 

specific heat at constant volume 

I 

specific heat at constant pressure 
constant (eq. (34)) 

function eissociated with second-order solution (eq. (28)) 
function used in description of mode (eqs. (lO)) 
function associated with second-order solution (eq. (3l)) 
steady rate of heat release per unit volume 
imaginary part 

Bessel function of first kind of order n. (ref. 6) 
constant of proportionality (eq. (Ub)) 
an5)litude (eq. (24)) 

axial extent of region of substantial temperature change 
(fig. 2) 

mean flow Mach number 

each of the sequence of numbers for which = 0 
mode 

Bessel function of second kind of order n (ref. 6) 
pressure coefficients (eqs. (37)) 

static pressure 
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ro.ri 

T 

t 

u 

X 

Z 

z 

a 

P 

r 

T 

e 

p 

T 

T* 

<P 

s 
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magnitude of disturbance velocity vector 

disturbance velocity vector 

inside radius of combiistion chamber (fig. l) 

radial coordinate (fig. l) 

combustor dimensions (fig. l) 

static ten 5 >erataire 

time 

axial perturbation velocity 

axial coordinate (sketch (a)) 

2R 

dimensionless quantity, — 

Of 

r 

dimensionless coordinate, — 

af 


2 ro 2r^ 
a£» ^ a^ 

coefficient (eq. (lOa)) 
coefficient (eq. (lOa)) 
ratio of specific heats, c^c^ 

an 5 )litude parameters (eqs. (2l) and (36), respectively) 
coefficient involved in solution of equation (l2) 
angular coordinate (fig. l) 
density 

time lag (eq. (l5)) 
most probable time lag (eq. (l7)) 
distvirbance velocity potential (q = Vcp) 
frequency in absence of burning 
frequency 
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Superscripts: 

(1) first-order solution 

(2) second-order solution 

' disturbance quantities and, -when convenient, ordinary 

differentiation 

Barred symbols denote flow quantities in absence of disturbance. 
Subscripts: 

c region rQ<r<r^ 

n mode 

w region r^<r<R 

Subscript notation for partial differentiation ^dien convenient. 
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Figure 1. - Aaaumed conTlguratlon of combustor. 
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Eadlal coordinate, r/B 

Figure 5. - Eadlal distritutlon of pressure amplitude. 
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Figure 6. - Modlfloatlon of fundamental pressure fluctuation at cylinder 
wall, owing to second-order terms. 
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